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1.  INTRODUCTION 


1.1  Background.  There  is  a  need  to  predict  blast  wave  overpressure  levels  accurately 
over  a  large  range  of  distances  from  the  gun  muzzle.  For  instance,  designing  an  enclosure 
for  reducing  impulsive  noise  requires  an  estimate  of  the  forces  and  impulses  on  its  inside 
surfaces.  The  enclosure’s  walls,  which  reduce  noise  emanating  from  large  guns,  may  be  as 
far  away  as  100  calibers  from  the  gun  muzzle.  But  it  is  aJso  important  to  know  muzzle-blast 
levels  closer  to  the  gun.  Muzzle  blast  can  injure  people  and  harm  material  at  shorter  ranges. 
It  is  therefore  important  to  be  able  to  predict  blast  levels  for  new  or  upgraded  weapons 
before  they  are  fired  to  avoid  unexpected  damage. 

The  muzzle  blast  levels  depend  on  the  weapon  design  and  emplacement,  the  propellant 
and  projectile  characteristics,  and  the  launch  conditions.  By  generalizing  and  extending 
Smith’s  work  (1974),  the  U.S.  Army  Research  Laboratory  (formerly  the  Ballistic  Research 
Laboratory  [BRL])  developed  prediction  methods  for  bare  muzzle  guns  based  on  data  col¬ 
lected  in  the  range  of  10  to  50  calibers  from  the  gun  muzzle  (Pansier  &  Schmidt  1983; 
Pansier  1985;  Heaps,  Pansier,  &  Schmidt  1985).  This  approach  was  also  based  on  the  cou¬ 
pling  between  solutions  to  the  problems  of  variable  energy  blast  waves  and  blast  waves  from 
asymmetrically  initiated  charges.  From  dimensional  analysis  (Baker  1973),  a  scaling  length, 

is  obtained.  This  scaling  length  depends  on  such  parameters  as  the  exit  muzzle  pressure, 
exit  temperature  of  the  propellcint,  and  the  field  angular  position  from  the  gun  bore  axis. 
The  field  distance  from  the  muzzle  divided  by  a  scaling  length  is  used  as  the  universally 
independent  variable  to  plot  the  peak  overpressure,  and  so  forth.  Free  parameters  are  gener¬ 
ated  by  the  dimensional  analysis  and  are  determined  by  a  least  squares  fit  to  the  data.  The 
resultant  predicted  free  field  muzzle  blast  is  implemented  on  a  computer  and  can  be  applied 
for  blast  waves  incident  on  surfaces  to  obtain  the  reflected  pressures. 

In  the  studies  cited  above,  data  were  collected  for  distances  close  to  the  muzzle.  The 
larger  distances  encountered  for  the  enclosure  generate  an  uncertainty  about  the  prediction 
accuracy  (especially  the  impulse).  Although  the  old  model  predicts  impulse  and  peak  over¬ 
pressures  close  to  the  muzzle,  an  examination  of  the  impulse  expression  shows  the  energy  of 
the  wave  is  increasing  although  it  should  decrease  with  larger  distances.  Other  people  have 
investigated  muzzle  blast  at  greater  distances.  Soo  Hoo  and  Moore  (1972)  primarily  studied 
various  naval  guns  with  data  taken  for  distances  between  20  and  110  calibers.  They  obtained 
pressure  contours  scaled  in  terms  of  calibers  for  naval  guns  with  bore  sizes  from  40  mm  to 
8  inch  (203  mm).  They  also  obtained  data  for  U.  S.  Army  20-mm  M3  and  M197  cannon. 
Naval  guns  of  different  bore  diameters  but  having  the  same  barrel  length  in  calibers  scaled 
well  over  a  large  range  of  distances.  Scaling  for  the  naval  guns  was  achieved  even  though 
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the  distance  in  calibers  from  the  ground  plane  to  the  muzzle  was  not  constant. 

Pater  (1981)  obtained  additional  data  in  the  far  field  for  naval  guns.  Using  Soo  Hoo  and 
Moore’s  data,  he  investigated  the  changing  relative  angular  distribution  of  the  gun  blast’s 
shock  wave  strength  with  distance.  The  peak  sound  pressure  level  (PSPL)  was  approximately 
23  dB  greater  at  the  front  than  at  the  rear  for  the  pressure  wave  front  located  20  calibers  from 
the  muzzle  of  a  5-inch/54  navzil  gun,  while  in  the  far  field,  the  PSPL  difference  decreased  to 
14.3  dB.  Schomer,  Little,  cind  Hunt  (1979)  used  guns  ranging  from  mortars  to  tank  cannon  to 
obtain  PSPL  differences  ranging  from  14  to  18  dB  in  the  far  field.  Most  recently,  Kietzman, 
Pansier,  and  Thompson  (1991)  have  obtained  data  from  a  105-mm  tank  cannon  and  noted 
that  the  angular  distribution  of  the  shock  wave  strength  changed  with  distance.  The  original 
model  of  Pansier  and  Schmidt  (1983)  assumed  an  unchanging  distribution  with  distance. 

1.2  Objectives.  It  is  the  purpose  of  this  investigation  to  develop  a  model  for  muzzle 
blast  that  accurately  predicts  quantities  for  a  lairge  range  of  distances  from  the  muzzle.  The 
physics  of  the  blast  wave  are  to  be  emphasized  while  constructing  the  model,  which  uses  a 
scaling  approach  with  free  parameters  to  be  fit  with  the  experimental  results.  To  these  ends, 
data  are  obtained  using  a  wide  variety  of  weapon  parameters  using  a  large  range  of  distances 
to  provide  an  adequate  data  base  for  the  improved  model. 

2.  EXPERIMENTAL  SETUP 

The  experiment  to  obtain  the  small  caliber  data  was  conducted  at  the  ARL  indoor 
aerodynamics  range  for  distances  15  to  100  calibers  from  the  muzzle.  Because  of  limited 
space  in  the  indoor  aerodynamics  range,  the  ARL  transonic  range  was  used  when  the  gauges 
were  placed  400  calibers  from  the  muzzle.  The  weapons  used  in  the  test  were  a  300-magnum 
barrel,  another  300-magnum  barrel  that  had  a  reduced  bore  length,  and  a  shortened  carbine 
barrel.  A  schematic  of  the  gauge  positions  around  the  gun  muzzles  is  shown  in  Pigure  1. 
In  addition,  a  gauge  was  positioned  near  the  muzzle  to  establish  a  zero  reference  time  for 
the  experiment.  The  current  small  caliber  data  were  supplemented  by  105-mm  tank  cannon 
data  obtained  in  the  Kietzman,  Pansier,  and  Thompson  (1991)  experiment. 

The  loadings,  velocity,  and  muzzle  pressure  axe  shown  in  Table  1.  The  muzzle  pressure 
in  the  leist  column  is  the  peak  value  just  before  the  projectile  exits  the  barrel.  With  the 
exception  of  the  last  row,  the  first  number  in  the  description  column  of  Table  1  refers  to 
the  projectile  mass  in  grains,  while  the  second  number  refers  to  the  charge  mass  in  grains. 
Purther  discussions  about  the  firings  will  use  these  descriptions  for  identification.  The  last 
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Figure  1.  Gauge  Positions  for  Study 

row  refers  to  the  parameters  for  the  105-mm  cannon  shooting  the  M735  round. 

3.  MODELING  APPROACH 

The  authors  predict  four  quantities  that  partially  characterize  the  blsist  wave.  Three  of 
these  quantities  are  illustrated  in  Figure  2.  The  peak  overpressure,  P  =  (pp  -  Poo  )/?<»,  is 
shown  as  is  the  time  of  arrival,  ta,  and  the  positive  phase  duration,  t.  The  time  of  arrival  is 
simply  the  time  elapsed  between  the  projectile  exiting  the  muzzle  and  the  arrival  of  the  wave 
at  the  selected  field  point.  The  positive  phase  duration  is  the  time  difference  between  the 
passage  of  the  wave  front  and  the  passage  of  the  wave’s  change  from  positive  overpressure  to 
negative  overpressure.  The  impulse  corresponding  to  the  time  integral  of  the  wave’s  positive 
phase  is  also  of  interest. 

For  gun  muzzle  blast,  the  energy  efflux  from  the  muzzle  declines  monotonically.  For 
instantaneous  energy  deposition,  dimensional  analysis  leads  to  scaling  relationships  that 
generate  universal  curves.  Gun  blast  can  be  treated  in  a  similar  manner  but  in  place  of  the 
instantaneous  energy,  E,  the  peak  energy  efflux,  dEldt,  is  used.  The  peak  overpressure, 
P  =  Pp  —  Poo,  is  expressed  in  functional  terms  as, 

P  =  P{r,p^,aoo,  dEldt),  (1) 
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Barrel 

Description 

Propellant 

Type 

Charge 
Mciss 
g  (gr) 

Projectile 
Velocity 
m/s  (ft/s) 

Muzzle 

Pressure 

MPa  (kpsi) 

300 

Long 

Magnum 

1220w70 

4831 

4.82  (70.0) 

899  (2950) 

59.40  (8.62) 

1220wl2 

Bullseye 

0.83  (12.0) 

282  (925) 

15.56  (2.26) 

1125w75 

4831 

5.16  (75.0) 

975  (3200) 

80.65  (11.69) 

1125wl2 

Bullseye 

0.83  (12.0) 

457  (1500) 

20.37  (2.95) 

300 

Short 

Magnum 

s220w36 

4227 

2.48  (36.0) 

594  (1950) 

74.9  (10.85) 

s220w9 

Bullseye 

0.62  (9.0) 

221  (725) 

17.7  (2.57) 

s125w342 

4227 

2.89  (42.0) 

792  (2600) 

89.60  (12.99) 

C2krbine 

cl3p6 

2400 

0.94  (13.6) 

518  (1700) 

6.59  (9.56) 

105  mm 

105mm 

M30 

5966  (-) 

1501  (4925) 

71.5  (10.37) 

in  which 

Pp  =  peak  value  of  the  pressure  at  a  given  field  point, 
r  =  distance  from  muzzle  to  field  point, 

Poo  =  2anbient  density, 

Ooo  =  ambient  speed  of  sound, 
dE/dt  —  energy  deposition  rate  into  atmosphere. 
Here  the  energy  deposition  rate  can  be  written  as. 


dE 

dt 


7»-l 


[1  + 


(7e  -  1) 
2 


(2) 


in  which  Ae  is  the  area  of  the  bore.  Mg  is  the  exit  Mach  number  of  the  propellant  flow 
immediately  after  the  projectile  exits  the  muzzle,  pe  is  the  peak  muzzle  pressure  divided  by 
the  atmospheric  pressure,  Ug  is  the  velocity  of  the  exiting  propellant,  and  7g  is  the  specific 
heat  ratio  for  the  exiting  propellant.  For  supersonic  flow  the  exit  conditions  would  be 
identical  to  the  conditions  just  before  projectile  releeise  from  the  muzzle.  For  subsonic  flow 
the  exit  conditions  are  different  and  are  discussed  by  Fansler  and  Schmidt  (1983). 


Using  the  Buckingham  Pi  theorem,  a  scaling  length  for  constant  energy  efflux  explosion 
is  obtained. 


^  ~  \J{dEfdt)l{p^aoo) 


(3) 
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Figure  2.  Blast  Wave  Values  to  be  Obtained 
and 

P  =  P/p^  =  P(r/£).  (4) 

Also  the  time  of  arrival,  and  the  positive  phase  duration,  r,  scale  as 

taaoo/i  =  fir/i),  (5) 

2md, 

TOoo/i  =  9{r/i).  (6) 

The  above  treatment,  which  assumes  spherically  symmetrical  conditions,  is  modified  for 
gun  blast,  where  the  peak  overpressure  varies  strongly  with  the  polar  angle,  0,  from  the 
boreline.  Following  Smith  (1974),  we  obtain  a  scaling  length,  P,  that  depends  on  the  polar 
angle, 

Pji  =  fi  cos  9  yj\  —  (P  sin^  9,  (7) 

in  which  p  is  the  momentum  index,  which  determines  the  ratio  of  the  peak  overpressure  (with 
distance  held  constant)  to  the  front  along  the  boreline  divided  by  the  peak  overpressure  to  the 
rear  along  the  boreline.  With  the  modification,  we  have  for  the  gun  blast  peak  overpressure, 
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in  which  f  =  rjP. 


P  =  P{f). 


(8) 


3.1  Peak  Overpressure.  Before,  we  assumed  a  one-term  equation  that  was  fitted 
with  the  data  to  obtain  the  value  of  /x  and  the  exponent  power  of  f.  That  is,  P(0‘’)/P(180°) 
did  not  vary  with  the  distance  from  the  muzzle.  But  here,  a  nonlinear  equation  is  assumed 
so  that  the  above  ratio  of  peak  overpressures  changes  with  distance  from  the  muzzle.  The 
predictive  expression  models  an  acoustic  wave  at  large  distamces  but  closer  in  the  peak  over¬ 
pressure  versus  distance  relationship  has  a  steeper  slope,  as  also  occurs  with  instantaneous 
explosions.  The  following  expression  for  the  peak  overpressure,  P,  is  assumed: 


P  = 


(9) 


Equation  (9)  is  matched  with  pressure  data  and,  as  before,  uses  Equation  (7)  to  vary  the 
strength  of  the  blast  with  the  polar  angle,  B. 

Originally,  we  used  Equation  (9)  with  the  data  obt^ned  in  the  far  field  to  obtain  the 
value  of  the  momentum  index.  The  dependence  of  the  blast  strength  with  angle  would  thus 
be  well  established  in  the  far  field.  The  values  of  A  and  B  could  then  be  determined  by  least 
squares  fitting  with  both  the  near  and  far  field  data.  Since  it  is  well  known  that  the  value 
of  p  in  the  far  field  depends  on  a  particular  weapon,  significant  parameters  were  sought  for 
the  blast  strength  dependence  on  angle. 

One  candidate  is  the  blow-down  pa”ai:  .:t“r,  6,  which  is  essentially  the  ratio  of  the  scales 
for  the  gun  tube  emptying  time  and  the  time  scale  for  the  muzzle  blast 


8  = 


Lapp 

Vj,i  ' 


(10) 


in  which  L  is  the  length  of  the  gun  barrel  in  calibers,  a^p  is  the  ambient  sound  speed,  and 
Vp  is  the  muzzle  velocity  of  the  projectile.  Pansier  and  Schmidt  (1983)  used  a  numerical 
code  developed  by  Erdos  and  DelGuidice  (1975)  to  establish  the  dependence  of  the  peak 
overpressure  upon  8.  Figure  3  shows  how  the  peek  '.vernressure  varies  with  barrel  length, 
which  is  equivalent  to  varying  8.  Guided  by  the  results  of  the  Erdos  code.  Pansier  (1985) 
found  that  the  measured  positive  phase  duration  and  impulse  descriptions  could  be  improved 
by  including  the  parameter,  ^  =  StjP,  into  the  least  squares  fitting.  The  parameter  /?  has 
the  same  meaning  as  8  except  now,  the  angular  position  of  interest  is  taken  into  account. 
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Figure  3.  Variation  of  Peak  Overpressure  with  Barrel  Length 

The  location  of  the  Mach  disc  is  also  used  as  a  parameter  to  lit  the  dependence  of  blast 
strength  on  the  angle  as  is  the  scale  length  ifD  \n  which  D  is  the  diameter  of  the  bore. 
The  authors  also  used  an  apparent  source  point  for  the  gun  blast,  which  is  forward  of  the 
muzzle  at  a  fixed  fraction  of  the  distance  from  the  muzzle  to  the  Mach  disc  for  an  equivalent 
steady  jet.  Shadowgraphs  of  the  blast  wave  show  an  almost  spherical  structure  whose  center 
is  forward  of  the  muzzle.  For  the  steady  jet,  the  position  of  the  Mach  disc  relative  to  the 
muzzle  for  the  steady  jet  is  given  as 

x'mID  =  A/e^7epe/2.  (11) 

For  this  study,  data  are  fitted  for  the  source  point  position,  Xo  =  O.STSxJi^ / D,  and  for  Xo  =  0. 

The  values  of  the  significant  parameters  are  given  in  Table  2.  The  dependence  of  the 
momentum  index,  fi,  is  explored  with  these  parameters. 

The  authors  first  used  the  approach  given  above  on  the  data  for  distances  >100  cal¬ 
ibers.  The  blow-down  parameter,  was  investigated  because  of  the  results  obtained  in 
Figure  3,  which  shows  a  different  overpressure  fall-off  for  different  values  of  barrel  length,  or 
equivalently  /3.  The  scaled  distance  can  be  written  in  terms  of  some  parameters  of  interest. 
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Table  2.  Significant  Blast  Parameters 


Designation 

tfD 

8 

cl3p6 

9.44 

1.85 

15.80 

s220w9 

4.95 

18.10 

7.07 

s220w36 

9.76 

3.75 

17.50 

sl25w42 

10.70 

2.39 

22.20 

1220wl2 

3.47 

41.30 

5.33 

1125wl2 

3.66 

23.00 

6.40 

1220w70 

8.73 

4.78 

22.20 

1125w75 

10.29 

3.60 

25.80 

105mm 

13.84 

0.85 

37.25 

rie  =  (TlD)(D/e)(tlOKII')-  (12) 

Here,  is  the  value  that  would  be  obtained  using  am  infinitely  long  tube  whereas  t'  is 
the  scale-length  value  that  depends  upon  the  field  position  for  tubes  of  finite  length.  The 
behavior  of  P  is  assumed  as 

^  =  ^{1  -  [1  -exp(-Bif)]exp(-yl,^)},  (13) 

in  which  =  81/ P.  For  large  distances,  the  scaled  length  would  asymptotically  approach 
the  ratio  of  P  jP^  corresponding  to  the  blow-down  parameter  value.  For  small  distances,  the 
differences  for  long  tubes  and  short  tubes  should  diminish  since  the  peak  values  would  occur 
before  the  gun  tube  had  emptied  appreciably.  The  exponential  term  involving  the  distance 
is  used  to  allow  for  this  expected  phenomenon. 

3.2  Time  of  Arrival.  As  in  Fansler  and  Schmidt  (1983),  the  pressure-jump  Mach 
relation. 


can  be  equated  to  the  predictive  equation  for  peeik  overpressure  and  integrated  to  obtain  a 
closed  form  expression  for  the  time  of  arrival, 
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h  -  h  =  A-Cf)  -  yin [2A'(f)  +  2f  +  A'], 


(15) 


in  which, 


t  = 


tCLoo 

I' 


X(r)  =  >/f2  +  A'f  +  B'  , 


A'  = 


27A 
7  +  1  ’ 


fl'  = 


2'iB 

7  +  1 


2uid  7  is  the  specific  heat  ratio  for  ambient  air.  The  quantity  A/,  is  the  Mach  number  for 
the  shock  front  moving  through  the  air. 


3.3  Positive  Phase  Duration.  The  wave  front  travels  at  the  speed  determined  by 
the  peak  overpressure  of  the  wave.  The  point  in  the  wave  where  the  overpressure  becomes 
zero  travels  at  the  speed  of  sound.  The  time  for  its  arrival,  ic-,  is 

fo  =  f  +  G.  (16) 

Here  G  is  a  value  that  depends  on  other  parameters,  for  instance,  /?.  The  positive  phase 
duration  is 


f  =  ic-h,  (17) 

or  substituting  in  the  expressions  for  and  ic-, 

f  =  f-  X(f)  +  y  In  [2X(f)  +  2f  +  at  -  +  G.  (18) 

The  positive  phase  duration  data  will  be  used  to  find  a  form  for  G  that  will  minimize  the 
root  mean  square  error. 

3.4  Impulse.  The  impulse  I  is, 

1=  f\p-P^)dl.  (19) 

The  dimensionless  form  of  the  impulse  is. 
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A  form  for  the  impulse  is  sought  that  is  consistent  with  the  present  treatments  of  the  peak 
overpressure,  time  of  arrival,  and  positive  phase  duration.  The  Friedlander  waveform  has 
been  used  successfully  to  describe  the  shape  of  the  initial  positive  phase  component  of  the 
wave  for  blasts  (Baker  1973).  The  most  primitive  form  is, 

(p-  /*oo)/Poo  =  P[l  -(t-  <a)/rlexp[-(<  -  t„)/r].  (21) 

Assuming  the  Friedlander  relationship  for  the  positive  part  of  the  waveform  in  Equation  (19), 
we  obtain  the  simple  expression, 

7  =  Pf/e.  (22) 

To  obtain  a  fit,  we  use  the  impulse  data  with  the  already  fitted  form  of  Equation  (9)  and 
with  Equation  (18),  with  G  to  be  determined  by  least  squares  fitting  of  the  impulse  data. 

4.  RESULTS 

4.1  Peak  Overpressure.  First,  a  fit  was  made  to  the  far  field  data  obtained  at  100 
and  400  calibers.  The  results  of  the  fits  with  various  parameters  used  are  given  in  Table  3 
in  which  A  is  the  multiplying  factor  in  Equation  (9).  At  these  distances,  the  contribution  of 
the  second  term  in  Equation  (9)  is  insignificant.  The  right-hand  column  gives  the  root  mean 
square  (rms)  error  for  the  fit. 

The  supposition  that  a  source  location  in  front  of  the  muzzle  improves  the  data  fit 
contradicts  our  findings.  However,  other  approaches  are  more  successful.  The  rms  error 
vzdue  is  reduced  when  the  momentum  index,  ft,  is  allowed  to  vary  with  Mach  disc  location. 
The  use  of  7/7?  as  a  parameter  also  provides  good  correlation.  That  both  parameters  provide 
good  correlation  is  not  surprising  since  the  two  parameters  increase  as  the  square  root  of  the 
muzzle  exit  pressure. 

The  question  whether  £'  depends  on  ^  was  investigated  by  fitting  Equation  (13)  to  the 
results  of  Model  10  in  Table  3.  The  terms  within  the  brackets  approach  unity  in  the  far  field 
and  were  set  to  unity  in  the  fitting  procedure.  The  parameter  Ai  is  an  insensitive  parameter 
for  fitting  and  so  the  blow-down  parameter  was  abandoned  as  a  parameter  to  be  used  for 
any  further  fitting  studies  of  the  peak  overpressure. 


Table  3.  Least  Squares  Fit  Results  for  Far  Field  Peak  Overpressure  Data 


Model 


Momentum  Index  (/x) 

HI 

Source  Position 

RMS 

0.71 

0.88 

0.0 

0.311 

0.69 

0.89 

0.375x'w/D 

0.324 

0.82 

0.93 

6(0.375x'jvf/D) 
b  =  -4.8 

0.271 

a  +  b{e/D) 

0  =  0.91,6=  -0.031 

0.90 

O.Z75x'm  ID 

0.270 

o  +  b{£/D) 
a  =  0.9, 6  =  —.026 

0.90 

0.0 

0.267 

o  +  6(£/D)  +  c6 
a  =  .88,6=  -.028,  c  =  .007 

0.90 

0.375x'jv,/D 

0.271 

o  +  b{i/D)  +  c6 
a  =  .87,6  =  -.023,c  =  .007 

0.90 

0.0 

0.268 

a-\-  bS 

a  =  .61,6=  .0054 

0.89 

0.375x',^/D 

0.281 

0  +  66 

a  =  .65, 6  =  .0046 

0.90 

0.0 

0.275 

o  +  6x'a//T> 
o  =  .87, 6  =  —.01 

0.91 

0.0 

0.263 

o  +  b{x'^ID)IM. 
a  =  .87,6=  -.012 

0.91 

0.0 

0.267 

1' ft  =  a6cos  0  -k-  y/\  —  a?  sin^  0 

0  =  0.31,6=2.4 

0.76 

0.0 

0.287 

A  similar  parameter  study  was  performed  with  all  the  data.  The  results  of  the  far  field 
study  were  used  to  fix  some  of  the  parameters  for  some  of  the  models.  With  these  parameters 
fixed,  the  accuracy  in  the  fM  field  is  emphasized  at  the  expense  of  the  predictions  for  the 
near  field.  Table  4  gives  the  value  of  together  with  the  values  of  the  multiplying  factors, 
A  and  B,  obtained  with  various  approaches. 

In  addition,  a  fit  was  made  to  the  data  as  in  the  previous  work  of  Fansler  and  Schmidt 
(1983).  With  this  set  of  data,  the  results  were, 

P  =  1.88{-J  ,  (23) 

in  which  fi  =  0.72.  This  compares  with  the  old  data  fit  of 
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Table  4.  Least  Squares  Fit  Results  for  All  PeeJc  Overpressure  Data 


Model 

Momentum  Index  (/i) 

tm 

B 

Source  Position 

RMS 

1 

0.7,  Fixed 

0.80 

1.99 

0.0 

0.283 

2 

0.71 

0.82 

iic?a 

0.375x'm/D 

0.283 

3 

a  +  b{(fD),  Fixed 
a  =  0.91, 6=  -0.031 

0.88 

1.50 

0.375x'mID 

0.354 

4 

a  +  b{eiD) 

0  =  0.90,6=  -0.026 

0.84 

1.55 

Q.375x*m/D 

0.267 

5 

o  +  b{t/D) 
a  =  0.81,6=  -.019 

0.89 

1.56 

0.0 

0.256 

6 

a  + 

a  =  .85,6  =  —.01 

0.83 

1.57 

0.375x'a,/£) 

0.260 

7 

«  +  b{x'iy(/D) 

a  =  .83,6=  -.0063 

0.89 

1.61 

0.0 

0.249 

8 

a  +  6(x^/Z?)/A/e,  Fixed 
a  =  M,b=  -.012 

0.84 

1.61 

0.0 

0.268 

9 

0  +  6(z'jv,/i))/A/e 
a  =  .85,6=  -.0086 

0.88 

1.59 

0.0 

0.259 

in  which  =  0.78.  The  value  of  the  momentum  index  is  smaller  for  the  new  data  and  reflects 
the  influence  of  the  data  in  the  far  field  in  which  the  pressure  differences  with  angle  become 


From  this  study,  the  predictive  fit  selected  was  model  7,  where  n  depends  upon  the  Mach 
disc  location.  With  A  and  B  determined,  the  working  version  of  Equation  (9)  can  be  written 


P  =  .89(^)+l,6l(^y,  (25) 

in  which 


H  =  0.83  -  0.0063(x;j^/£>).  (26) 

The  result  for  (i  shows  that  the  ratio  of  the  front  to  back  shock  strengths  becomes  lower 
with  higher  muzzle  exit  pressure.  This  result  agrees  with  the  data  of  Schomer  and  Raspet 
(1984).  Predictions  from  Equation  (25)  are  compared  with  the  correlating  data,  as  shown  in 
Figure  4. 
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SCALED  DISTANCE 

Figure  4,  Pressure  vs.  Scaled  Distance  with  New  Approach 

4.2  Time  of  Arrival.  The  time-of-aa'rival  expression,  Equation  (15),  can  be  used  with 
the  determined  constants  A  and  B  and  the  time-of-arrival  data  to  obtain, 


=  X(f)  -  0.52  In  [2X{f)  +  2f  +  1.04]  -  0.56, 


(27) 


in  which. 


A'(f)  =  Vf*  +  1.04f  +  1.88  . 


The  data  were  logarithmically  weighted  to  minimize  the  relative  error.  The  constants  ob¬ 
tained  with  no  weighting  have  comparable  values.  Least  squares  fits  were  also  tried  by 
re-determining  the  parameters  previously  found  with  the  peak  overpressure  data.  The  ob¬ 
tained  rms  error  was  not  improved  significantly  with  the  time-of-arrival  data  determining 
A,  B,  and  ft. 


Figure  5  shows  the  calculated  minus  the  observed  (residuals)  as  a  function  of  the 
scaled  length.  Approximately  700  data  points  are  shown  in  this  graph.  Most  of  the  points  are 
bunched,  with  a  relatively  few  outlying  data  points  scattered  about.  The  time-of-arrival  data 
were  obtained  with  the  aid  of  a  code  that  determines  automatically  the  peak  overpressure, 
time  of  arrival,  positive  phcise  duration,  and  impulse.  It  is  possible  that  some  of  the  outlying 
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SCALED  DISTANCE 

Figure  5.  Calculated  h  Minus  the  Observed  fa  (residuals)  vs.  r/£' 

data  with  negative  values  could  be  caused  by  triggering  on  the  precursor.  The  origins  of 
errors  in  the  data  have  not  yet  been  thoroughly  investigated.  Nevertheless,  most  of  the  data 
aire  densely  located  near  the  zero  line  except  for  the  larger  scaled  distances.  Another  source 
of  error  is  the  uncertainty  of  the  sound  speed.  The  experiment  was  conducted  covering  a 
period  from  April  to  September,  1992,  and  significant  differences  in  temperature  (and  thus 
sound  speeds)  could  have  occurred  over  the  times  when  data  were  taken. 

4.3  Positive  Phase  Duration.  The  positive  phase  duration  data  were  obtained  by 
using  the  code  mentioned  above.  The  code  tests  for  the  sign  cf  the  value  of  overpressure 
to  change  from  positive  to  negative.  This  change  in  sign  is  also  related  to  obtaining  the 
impulse.  If  the  sign  changes  and  again  goes  positive  within  a  given  time  interval  and  a 
significant  amount  of  impulse  is  added,  another  sign  change  is  needed  to  indicate  that  the 
positive  phase  of  the  wave  has  been  completed.  Various  forms  for  G  were  assumed  to  find  a 
minimum  value  for  the  rms  error.  Table  5  gives  the  forms  tried  and  the  results. 


Table  5.  Positive  Phase  Duration  Parameter  Fitting  Results 


Model 

G 

6, 

62 

h 

EH 

RMS 

1 

6i 

1.22 

0.777 

warn 

61—62  cos  9 

0.572 

1.01 

0.636 

(■Ki 

61  +  63(1  —  exp(— 62/^)]  —  64  cos  ^ 

-  0.157 

.00264 

16.8 

8.14 

0.620 

l■■l 

61  +  62/?  —  63  cos  9 

0.0944 

0.0319 

0.584 

1 

0.583 
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4.4  Impulse.  As  mentioned  when  discussing  the  positive  phase  duration  results,  a 
computer  code  developed  to  reduce  the  data  has  the  capability  to  integrate  the  pressure 
data  to  obtain  the  impulse.  It  can  also  calculate  the  impulse  for  double  peaks  in  quick 
succession,  even  though  there  might  be  a  small  interval  in  which  the  wave  goes  negative. 
For  calculations  of  structurad  vulnerability  the  impulse  is  the  important  quantity.  For  our 
approach  discussed  earlier,  f  has  to  be  determined  in  Equation  (22).  Equivalently  for  this 
approach,  the  value  of  G  is  determined. 

The  seeded  impulse  as  a  function  of  scaled  distance  is  shown  in  Figure  6.  The  scatter 
is  more  pronounced  than  for  the  scaled  peak  overpressure.  The  two  lower  points  were  later 
eliminated  as  being  processed  improperly  by  the  computer  code. 


SCALED  DISTANCE 

Figure  6.  Scaled  Impulse  vs.  Scaled  Distance 

Table  6  gives  some  of  the  models  tried  in  terms  of  functions  of  G  and  the  results.  Model  19 
is  chosen  to  calculate  the  positive  phaise  duration  and,  using  Equation  (25),  the  impulse.  The 
value  of  G  is, 

G  =  0.09  -  0.003796  +  1.07[1  -  1.36  exp  (-0.049f)]£/6'.  (28) 

Figure  7  shows  the  scaled  impulse  plotted  with  the  predicted  values  for  the  firings  that 
used  the  220  grain  projectile  with  12  grains  of  powder  for  the  long  barrel  (designated  as 
1220wl2).  Figure  8  gives  corresponding  results  for  the  sl25w42  firings.  These  figures  provide 
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Model 


1 


2 


3 


4 


Table  6.  Impulse  Parameter  Results 


62  f>3 


61  + 


Oi  +  62 

+b30 


61  +  63/5 
— 64003^ 


61—62  cos  6 


61  +  62^3 
—63  cos  d 


61  +  62O  —  63  cos  a 


61  +  63(1  -  exp  (-62/?)] 
—64  cos  6 


61  +  hP 

—63  cos  0  —  64/3  cos  0 


61  +  63 


61  +  62^6+63^/^' 


6i  +  62^/8  + 


61  +  626  +  bsi/c 
61  +  62^3  +  63^/ 


61  +  62(1  —  exp(— 63f) 
(1  -  €Xp(-b40) 


0.270 


0.0067 


0.352 


-.0649  0.0307 


-.509  0.539 


-.488  0.261 


0.134  i  0.972  I  0.0272  I  .455 


0.556  0.842 


0.146  0.024  0.459 


0.310  1  0.0276  I  0.744 


0.120 


9.39 

0.446 

0.411 

0.223 

-.137  0.00781  5.78 


0.608 


-.904  0.162 


-1.75  0.174 


-.680  0.193 


-.486  0.0189 


-.100  4.179 


1.49 


0.619 


0.413 


0.0868  0.0179 


-.268  0.0159  0.136  29.1 


61  +  62/3 


61  “f"  630 

+63(1  -  h exp {-b^f )]£/£') 


-.0085 


0.0337 


0.108 


1.07 


0.049 


RMS 


0.538 


0.440 


0.435 


0.435 


0.418 


0.468 


0.418 


0.450 


0.417 


0.438 


0.421 


0.428 


0.425 


0.410 


0.424 


0.433 


0.358 


results  from  contrasting  firing  conditions.  The  1220wl2  firings  constitute  firing  from  a  long 
barrel  with  low  velocities;  the  sl25w42  firings  constitute  firing  from  a  short  barrel  with  high 
velocities.  At  the  greater  distances,  the  scaled  impulse  values  are  higher  for  the  larger  polar 
angles.  Just  the  opposite  behavior  occurs  at  the  smallest  distances  and  the  predictive  curves 
reflect  these  trends  also.  We  cannot,  with  certainty,  explain  these  data  trends.  However,  we 
speculate  that  the  frontal  part  of  the  wave  may  be  feeding  energy  around  the  positive  phase 
of  the  wave  to  the  rear  of  the  gun. 

A  more  careful  analysis  should  also  consider  how  the  wave  shape  changes  with  distance 
and  angle.  Wilcoski  and  Pater  (1990)  show  that  the  relative  spectral  frequency  content  of 
the  wave  changes  significantly  with  distance.  One  could  take  this  change  into  zwrcount  with 
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Figure  7.  Comparison  of  Scaled  Impulse  with  Predicted  Values  for  1220wl2  Firings 


the  use  of  the  the  modified  Friedlander  wave, 


(p  -  Foo)/poo  =  P[1  -  (f  -  ia)lr]  exp  [-hit  -  fa)/r].  (29) 

Here,  h  is  the  wave  shape  factor  that  varies  with  distance. 

Data  for  instauitaneous  explosions  (Baker  1973)  show  that  the  value  of  h  declines  with 
distance  to  a  local  minimum  with  a  value  near  2.  Since  the  impulse  is  used  more  than  the 
positive  phase  duration  for  structural  computations,  we  have  concentrated  on  getting  a  good 
fit  to  the  impulse.  For  a  constant  value  of  impulse  and  peak  overpressure,  the  normalized 
positive  phase  duration  as  a  function  of  the  wave  shape  factor,  h,  is  shown  in  Figure  9. 
The  positive  phase  duration  is  normalized  by  its  value  when  h  =  1.  For  our  purposes,  the 
assumption  of  h  =  1  will  only  adversely  affect  the  time  representation  of  the  directed  wave 
and  reflected  waves  that  are  close  to  each  other. 

4.5  Predictive  Equations  for  the  Far  Field.  If  the  peak  overpressure  described  by 
Equation  (25)  is  applied  with  the  expression  for  the  positive  phase  duration  to  the  far-field, 
the  symbolic  and  mathematical  analysis  package.  Wolfram  Research’s  Mathematica  (1991), 
can  be  used  to  show  that  the  energy  in  the  blast  wave  has  a  logarithmic  singularity  at 
infinity.  Clearly  this  result  does  not  correspond  to  reality.  Landau  and  Lifshitz  (1959)  show 
dissipation  at  large  distances  from  the  explosion  results  in  the  peak  overpressure  declining 
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as, 


P  =  A2l[fyJ\n{fla2)].  (30) 

Using  Equation  (30),  one  can  obtain  dijdf  as  before  and  expand  the  result  in  series  form 
for  large  f.  Integration  to  obtain  the  time  of  arrival  and  then  the  positive  phase  duration 
shows  that  the  positive  phase  duration  will  grow  as  the  square  root  of  the  log  of  the  distance. 
The  energy  of  the  wave  will  then  decline  as  the  inverse  square  root  of  the  dist«ince.  It  was 
decided  to  match  Equation  (25)  and  Equation  (30)  at  f  =  50.  At  this  distance,  the  ratio  of 
blast  strengths  for  different  angul^  positions  is  changing  very  slowly  with  distance  and  the 
wave  behavior  is  almost  linear  as  for  an  Jicoustical  wave.  With  the  slopes  of  the  equations 
set  equal  to  each  other  at  this  distance,  the  constants  in  Equation  (30)  give  the  particular 
working  asymptotic  form  for  the  peak  overpressure. 


P  =  3.48975/[fy'ln(33119f)]. 


(31) 


In  turn  the  asymptotic  form  for  the  positive  phase  duration  is  also  found  and  is, 


f  =  2.99y'ln(33119f)  -  8.534  +  G. 


(32) 
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Figure  9.  Normalized  Positive  Phase  Duration  vs.  A,  a  Wave  Shape  Factor 

Figure  10  shows  a  comparison  of  the  asymptotic  form  for  the  duration,  Equation  (32), 
with  the  form  from  the  least  squares  fit.  Although  the  differences  are  small,  only  the  asymp¬ 
totic  forms  give  the  correct  energy  results  in  the  limit  of  large  distances. 

With  the  complete  solution  now  obtained,  we  can  represent  the  positive  phase  duration 
as  a  function  of  scaled  distance  and  angle,  which  is  shown  in  Figure  1 1  .  Although  increasing 
comparatively  rapidly  at  the  smaller  distances,  the  positive  phase  duration  increases  only 
slowly  with  larger  distances.  If  data  were  taken  only  close  to  the  weapon  for  a  limited  range 
of  r,  the  scaled  duration  as  a  functioii  of  r  would  have  its  smallest  values  at  small  f  and 
and  the  largest  values  of  r  would  be  fo  uid  at  a  larger  f  and  9.  The  equivalent  path  over  the 
surface  in  Figure  11  might  approximate  an  almost  linear  relationship.  Indeed,  Fansler  and 
Schmidt  (1983)  obtained  such  a  linear  fit. 

4.6  Predictive  Relationship  for  the  Complete  Wave.  When  waves  are  reflected 
from  a  surface,  the  direct  wave  and  reflected  wave  will  interact  at  some  points  in  the  field. 
The  front  of  the  reflected  wave  could  also  interact  with  the  negative  ph2ise  of  the  direct  wave. 
To  accurately  predict  the  wave  interaction,  a  knowledge  of  the  negative  phase  is  also  needed. 
Formerly,  the  Reed  wave  description  (1977)  was  used  for  both  the  positive  and  negative 
phase,  but  we  have  used  the  Friedlander  description  for  the  positive  phase.  The  Friedlander 
wave  is  only  descriptive  of  the  positive  phaae.  Thus,  we  need  to  join  the  Friedlander  positive 
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Figure  10.  Parametric  and  Asymptotic  Forms  for  the  Positive  Phase  Duration 


phase  to  the  negative  phase  of  the  Reed  description.  The  Reed  wave  is 


(P  -  =  p\\  -  {(  -  <.)/rl(l  -  ((  -  (,)/(3.67r)J(l  -  |((  -  (.)/(3.67t))^}.  (33) 

For  t  —  ta>  T,  this  equation  multiplied  by  a  factor  to  make  the  integral  of  the  complete 
wave  equal  to  zero  describes  the  negative  phase.  If  the  time  integral  for  the  complete  wave 
were  not  equal  to  zero,  the  total  momentum  of  the  wave  would  increase  with  distance  with 
no  limiting  bound.  The  required  factor  is  equal  to  0.751.  Figure  12  is  a  plot  of  the  resulting 
wave.  A  discontinuity  occurs  at  the  juncture  of  the  positive  phase  and  the  negative  phase. 
The  disparity  in  slope  would  be  even  greater  if  the  shape  factor,  h,  for  the  Friedlander  wave 
were  set  equal  to  two,  as  one  might  prefer  to  optimadly  simulate  the  wave  at  larger  distances. 

4.7  Summary  of  Working  Equations.  These  working  equations  are  used  to  im¬ 
prove  the  present  computer  implementation  for  the  muzzle  blast  (Heaps,  Fansler,  and  Schmidt 
1985).  The  peak  overpressure  for  f  <  50  is, 

P=,89(g+1.6l(^y, 

in  which  the  scaled  length  taking  the  polar  angle  from  the  boreline  into  account  is, 
f'/f  =  fi  cos  $  d-  ^1  —  /i*  sin*  9  , 
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Figure  11.  Scaled  Positive  Phase  Duration  vs.  Scaled  Distance  and  Angle 
and  n,  the  momentum  index  is, 

/i  =  0.83-0.0063(z'A//^)  • 

The  peak  overpressure  for  f  >  50  is, 

P  =  3.48975/[fyin(33119f)]. 

The  time  of  arrival  is, 

=  A’(f )  -  0.52  In  [2X(f)  +  2f  +  1.04]  -  0.56  , 
in  which, 

X{f)  =  >/f2  +  1.04f  +  1.88  . 

The  positive  phase  duration  for  f  <  50  is, 

f  =  f  -  A’(f)  +  0.52  In  [2A'(f)  +  2f  +  1.04]  +  0.56  +  G  , 
in  which, 

G  =  0.09  -  0.003796  +  1.07[1  -  1.36exp(-0.049r)]£/r  . 

The  positive  phase  duration  for  r  >  50  is, 
f  =  2.99yin(33119f)  -  8.534  +  G  . 


The  expression  for  the  impulse  is. 
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Figure  12.  Pressure-Time  Relationship  for  the  Complete  Wave 


/  =  Pf/e  , 

and  the  impulse  can  be  determined  from  the  above  working  equations.  The  overpressure-time 
relationship  over  the  positive  phase  is, 

(p  -  i’oo)/?*,  =  P[1  -(t~  ta)/r]  exp  [-(t  -  ta)/r]  . 

The  overpressure-time  relationship  over  the  negative  phase  is, 

(p  -  i’oo)/Pco  =  0.751P[1  -  (t  -  ta)/r][l  -(t-  t,)/(3.67r)]{l  -  [{t  -  t,)/(3.67T)]2}  . 

5.  SUMMARY  AND  CONCLUSIONS 


Overpressure  data  were  obtained  for  a  wide  range  of  distances  and  cingles  from  the  muzzle. 
The  data  were  then  used  in  a  parametric  investigation  to  improve  a  predictive  method 
based  on  dimensional  analysis.  In  addition  to  the  fundamental  length  scaling  parameter, 
other  parameters  investigated  were  the  Mach  disc  location  and  the  blow-down  parameter. 
Investigations  in  the  far  field  showed  that  the  blow- down  parameter  was,  for  the  range 
of  weapon  parameters  investigated,  not  a  significant  parameter  for  use  in  describing  peak 
overpressure.  Although  the  Mach  disc  location  could  not  be  used  to  improve  the  location  of 
the  source  point,  it  was  a  meaningful  parameter  for  helping  to  characterize  the  momentum 
index,  fi.  A  two-term  overpressure  equation  was  zissumed  to  take  into  account  the  nonlinear 
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blast  wave  near  the  muzzle.  The  new  peak  overpressure  equation  was  used  together  with  a 
result  from  gas  dynamics  to  obtain  a  closed  form  expression  for  the  time  of  arrival.  Only 
one  additional  parameter  needed  to  be  determined  by  legist  squares  fitting  to  obtain  the 
predictive  equation  for  time  of  arrival,  whereas  in  the  older  model,  two  parameters  were 
needed.  The  positive  phase  duration  is  constructed  by  subtracting  the  wave’s  time  of  arrival 
from  the  time  for  the  wave’s  zero  pressure  point  to  arrive,  assuming  that  this  point  travels 
at  the  speed  of  sound.  A  Friedlander  wave  is  assumed  with  a  value  of  the  wave  shape  factor 
equal  to  one.  With  the  assumed  shape,  the  Friedlander  wave  is  integrated  to  obtain  the 
impulse.  The  impulse  depends  only  upon  the  positive  phase  duration  and  the  peak  pressure. 
The  free  parameters  assumed  for  the  positive  phase  duration  are  then  fitted  to  the  impulse 
data. 

For  longer  distances,  a  smooth  join  was  made  to  the  fitted  pressure-distance  curve  with 
an  inverse  distamce  divided  by  the  square  root  of  the  logarithm  of  the  distance.  The  positive 
phase  duration  is  also  modified  to  grow  with  the  square  root  of  the  logarithm  of  the  distance. 
With  these  modifications,  the  total  energy  of  the  wave  decays  at  long  distances  as  the  inverse 
square  root  of  the  logarithm  of  the  distance.  Although  a  physically  valid  solution  is  obtained 
at  longer  distances,  use  of  this  prediction  method  should  be  limited  to  distances  where  the 
wave  does  not  extensively  interact  with  noise-absorbing  surfaces  or  is  not  subject  to  wave 
distortion  by  heterogeneous  atmospheric  conditions. 

For  detailed  study  of  the  interaction  of  the  direct  wave  and  the  reflected  wave,  the  negative 
phase  of  the  wave  is  also  needed.  The  Friedlander  wave  can  give  a  good  approximation  over 
the  positive  phase  of  the  wave  but  cannot  describe  the  negative  phase.  To  complete  the  wave, 
a  modified  description  of  the  negative  phase  of  a  Reed  wave  is  joined  to  a  Friedlander  wave. 
The  negative  phase  of  the  wave  conforms  with  the  necessary  condition  that  the  time  integral 
of  a  spherical  wave  equals  zero,  else  the  energy  of  the  wave  would  grow  with  increasing 
distance. 

Further  work  should  include  finding  the  wave  shape  factor  as  a  function  of  scaled  distance. 
The  value  of  the  impulse  is  not  affected  by  errors  in  representing  the  wave  shape  but  better 
prediction  of  wave  interactions  depends  upon  wave  shape  accuracy. 
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Research  Triangle  Park,  NC  27709-2211 

1  Commander 

Aviation  Applied  Technical  Dir. 

ATTN:  SAVRT-TY-MSMA,  G.  Moffatt 
Fort  Eustis,  VA  23604-5577 

3  Department  of  the  Army 

Construction  Engineering 
Research  Laboratory 
ATTN:  CERL-SOI, 

P.  Schomer 
L.  Pater 
J.  Wilcoski 
P.  0.  Box  4000 
Champaign,  IL  61820 

1  Commander 

ASD/YHT 

ATTN:  LTJ.  Palumbo 
D.  Curley 

Eglin  AFB,  FL  32542 

1  Commander  (Code  3433) 

Naval  Warfare  Center 
ATTN:  Tech  Lib 
China  Lake,  CA  93555 
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2  Commander 

Naval  Surface  Warfare  Center 
ATTN;  6X, 

J.  Yagla 
G.  Moore 
Dahlgren,  VA  22448 

1  Commander  (Code  6120C) 

Naval  Ordnance  Station 
ATTN:  Susan  Peters 
Indian  Head,  MD  20640 

1  Commander 

US  Naval  Air  Systems  Command 
ATTN:  AIR-604 
Washington,  DC  20360 

1  Commander  (Code  3892) 

Naval  Warfare  Center 
ATTN:  K.  Schadow 
China  Lake,  CA  93555 

1  Commander  (Code  730) 

Naval  Surface  Weapons  Center 
Silver  Spring,  MD  20910 

1  Director 

NASA  Scientific  &  Technical 
Information  Facility 
ATTN:  SAK/DL 
P.  O.  Box  8757 
Baltimore/ Washington 
International  Airport,  MD  21240 

1  McDonnell  Douglas 
ATTN;  Joseph  Smuckler 
1014  Ferngate  Lane 
Creve  Coeur,  MO  63141 

2  Alliant  Techsystems  Inc. 

ATTN:  MS  MN  50-2060, 

T.  Melanger 
S.  Langley 

600  Second  Street,  Northeast 
Hopkins,  MN  55343 

1  S  4:  D  Dynamics,  Inc. 

ATTN;  R.  Becker 
7208  Montrico  Dr 
Boca  Raton,  FL  33433-6930 


2  AAI  Corporation 

ATTN:  T.  Stcisney 
J.  Hebert 

P.  O.  Box  126,  MS  100-405 
Hunt  Valley,  MD  21030-0126 

2  Aerojet  General  Corporation 

ATTN:  W.  Wolterman 

A.  Flatau 
P.  0.  Box  296 
Azusa,  CA  91702 

1  Lockheed  Aircraft,  Inc. 

ATTN;  J.  Brown 
J.  Perez 

P.  O.  Box  33,  Dept.  1-330/UPLAND 
Ontario,  CA  91761 

1  General  Electric  Armament 

&  Electric  Systems 
ATTN.  R.  Whyte 
Lakeside  Avenue 
Burlington,  VT  05401 

1  Franklin  Institute 

ATTN;  Tech  Library 
Race  &  20th  Streets 
Philadelphia,  PA  19103 

1  Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Johns  Hopkins  Road 
Laurel,  MD  20707 

2  Loral  Corporation 
ATTN:  S.  Schmotolocha 

B.  Axely 

300  N.  Halstead  St. 

P.  O.  Box  7101 
Pasadena,  CA  91109 

2  McDonnell  Dougltts  Helicopter  Co. 

ATTN:  D.  Van  Osteen 
R.  Waterfield 
Mail  Station  D216 
500  E.  McDowell  Rd. 

Mesa,  AZ  85205 
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1  FN  Manufacturing,  Inc. 
ATTN;  George  Kontis 
Post  Office  Box  24257 
Columbia,  SC  29224 

1  Red  Eye  Arms,  Inc. 

ATTN:  David  Byron 
Gunn  Station 
507  N.  New  York  Ave. 
Winter  Park,  FL  32789 

1  Old  Dominion  University 

Mathematics  Department 
ATTN:  Dr.  Charlie  Cooke 
Norfolk,  VA  23508 

1  Scitec  Inc. 

ATTN:  Alex  Zislan 
100  Wall  Street 
Princeton,  NJ  08540 

1  Los  Alamos  National  Lab. 

ATTN:  Thomas  Davis 
Group  WX-4 
MS  G787 

Los  Alamos,  NM  87545 

1  Olin  Corporation 

ATTN :  Stephan  Faintich 

PO  Box  222 

St.  Marks,  FL  32355 

1  Atlantic  Research  Corp. 

ATTN:  Mark  Friedlander 
5945  Wellington  Road 
MS  G787 

Gainesville,  VA  22065 


Aberdeen  Proving  Ground 

2  Dir,  USAMSAA 

ATTN:  AMXSY-D, 

Mr.  W.  Brooks 
Mr.  R.  Conroy 

7  Cdr,  USACSTA 

ATTN:  STECS-TD,  Mr.  Kelton 
STECS-AAL,  M.  Maule 
STECS-AA-LA,  J.  Scheuren 
STECS-AS-LA,  S.  Walton 
STECS-DA,  P.  Paules 
STECS-RM-PF,  S.  Hinte 
STECS-AS-HP,  J.  Andrews 

2  Cdr,  USATECOM 

ATTN:  AMSTE-TE-R,  Mr.  Keele 
AMSTE-TA-R,  W.  Marshall 


1  Dir,  AMC 

Int.  Mat.  Eval.  Div. 

ATTN:  AMCICP-IM,  R.  Bloom 

2  Dir,  USARL 

ATTN:  AMSRL-HR-SD, 

G.  Garinther 
J.  Kalb 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1 .  ARL  Report  Number  ARL-TR-227 _ Date  of  Report  Septeri>er  1993 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  infonnation  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


DEPARTME^fT  OF  THE  ARMY 


OFROAl.  BUSINESS 


BUSINESS  REPLY  »L4IL 

FKT  CUSS  IliIMT  Ik  OOm.  WE.  MD 


Postage  will  be  paid  by  addressee 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 
Aberdeen  Proving  Ground,  MD  21005-5066 


